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Abstract 

The properties of amplitude modulated broadband Alfven waves is investigated. In particular, 
the dynamics of circularly polarized dispersive Alfven waves, governed by a derivative nonlinear 
Schrodinger equation, is analyzed using the Wigner formalism. The modulational instability of 
random phase dispersive pump Alfven waves is investigated, and it is shown that the spectral 
broadening gives rise to a new mode structure. 

PACS numbers: 52.35.Mw, 52.40.Db, 94.30.Tz 
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I. INTRODUCTION 

About thirty years ago, Rogister \1] introduced an elegant perturbation technique of 
the Vlasov-Maxwell system of equations in order to investigate amplitude modulation of 
magnetic field-aligned circularly polarized dispersive Alfven waves in an electron-ion plasma. 
He denved a derivative nonpar Schrodinger (DNLS) e q uario„, which is exacriy integrable 
by the inverse scattering transform method [2] . The DNLS equation for nonlinear dispersive 
Alfven waves have also been derived by means of the two-fluid - Maxwell system of equations 
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and the modulational instability of a constant a amplitude Alfvenpump 
wave has been investigated . Also reported is the algebraic envelope Alfven soliton [3|, |9j. 
The DNLS equation has been extended to include the effects of collisions and wave- 
10, The resulting modified DNLS equation has been used to 



particle interactions 



investigate damping of the envelope Alfvenic soliton. 

In this paper, we report on the nonlinear stability of broadband dispersive Alfven waves in 
plasmas. We introduce a Wigner-Moyal formalism for the modified DNLS equation 

and derive a kinetic equation for Alfvenic quasiparticles. While the DNLS equation is 
appropriate for studies of modulational instability of a coherent Alfven pump wave, the 
kinetic equation can be adopted to investigate the effects of partial coherence and random 
phases associated with broadband dispersive Alfven wave packets. 



II. GOVERNING EQUATIONS 

Let us consider the nonlinear propagation of circularly polarized dispersive Alfven waves 
along the external magnetic field B^z, where Bo is the magnetic field strength and z is the 
unit vector along the z-axis. The dynamics of modulated Alfven wave envelopes is governed 
by the modified DNLS equation |l0, [TJj 

id t u + d 2 z u + isd z (\u\ 2 u) = p[u, «*], (1) 

where the subscripts t and z denotes the corresponding derivatives. Equation (JTJ) governs 
the amplitude modulation of circularly polarized Alfven waves. Here u = (B x + iB y ) /(2\1 — 
/3| 1//2 .Bo), P = V 2 /Vjl, s = sign(l — /3), time and space coordinates are normalized by the 
ion gyroperiod u' 1 and the ion skin depth Va/uci — c/oj pi . For our purposes, we have in 
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the collision-dominated case 

P = iDd 2 z u, (2) 

where D = (l/2)[(rji/ p ) + c 2 7] 2 /4:7r]uj ci /V2, % is the coefficient of the ion viscosity and r} 2 is 
the resistivity. On the other hand, in a collisionless case, we have 

z' — z 

where P denotes the principal value. The coefficient C depends on the velocity distribution 
of the particle species. For Va S> V s , we have 

V. 



P 



c c 

i — d z (uH(z, t)) = i — d z I u P 



dz' 



(3) 
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<1, 



27rmi \ Va , 
where Vr e is the electron thermal speed. 

With (J2J and (jHJ), Eq. ((TJ) is a special case of the equation 

id t u + g d 2 z u + F(\u\ 2 ;t, z)u + id z [G(\u\ 2 ;t, z)u] = 0, 



(4) 



(5) 



where g is a complex valued constant and F and G are complex valued functions. This 
equation can be found in a wide variety of applications, such as Bose-Einstein condensation 
and nonlinear optics, apart from in plasma physics. Next we proceed to analyze Eq. (jHJ) for 
broad band waves. 

In order to study the stability of broad band waves governed by Eq. (JHJ), it is convenient 
to introduce the Wigner function 
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pit, z,k) = —J d( e ik t(u*(t, z + (/2)u(t, z - C/2)> 



(6) 



for the Alfven waves. Here the angular brackets represents the ensemble average. The 
Wigner function is given by the Fourier transform of the two-point correlation function, and 
as a generalized distribution function for quasi-particles it can be used to describe a broad 
band spectrum of the field represented by u, e.g. Alfven waves. We note that from © the 
equality 

/oo 
dkp(t,z,k) (7) 
-oo 

holds. Applying the time derivative to the Wigner function (jHJ) and using (jSJ), one obtains 
the kinetic equation 



d t p + 2 Re{g)kd z p + lm{g) {\d 2 z - 2k 2 ) p + 2 Im Fexp ( | d z d, 



P 



+d z \ Re 



Gexp | d z d k 
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Gexp ( | d z d k 



p = 0. 



(8) 



As an example of the versatility of Eq. (jHJ), we look at the case g = 1, F = I, and G = 0. 
Then Eq. (J3J) reduces to the regular NLSE, and Eq. (jHJ) becomes 



d tP + 2kd zP + 21 sin (j £S) p = 0, 



(9) 



an equation which to lowest order gives the Vlasov equation 



d t p + 2kd z p+{dj){d k p) = 0. 



(10) 



For Eq. (JTJ) we have F = 0, g = 1 — iD, and G = si — [C/Att)H, and the kinetic equation 
) takes the form 



d t p + 2kd z p — 2sA;i sin f | d z d k J p + s<9 z I cos d z d^j p = Lp, 
where the operator expression on the right hand side is defined according to 



Lp = D (\d 2 z - 2k 2 ) p 



in the collisional case (J2J), and 

C 

Lp = k 



in the collisionless case (J3J). 



lm(H) cos d z d k j + Re(H) sin ^| <9 2 <9 fc 
- Im(ff) sin (| d z (9 fe ) p + Re (if) cos <9 2 <9 fc ) p 



(11) 



(12a) 



(12b) 



III. ANALYSIS OF THE COHERENT CASE 



We note that Eq. has solutions of the modified plane wave form u(t, z) = 
u(t) exp[ik z — iuo(t)t]. In the collisional case (J2J) the time-dependent frequency is given 
by 

Mt) = k 2 Q + - exp(-2Dk 2 t)}, (13a) 

while the amplitude is exponentially damped according to 



u{t) = u exp(—Dklt). 



(13b) 



Here Uq is the constant amplitude of the solution, and we note that as D — > we obtain the 
time-independent dispersion relation u = k^ + sk ul- In the collisionless case (J3J) however, 
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the character of the function p[u, u*] together with the above plane wave ansatz makes the 
amplitude time-independent, and subsequently the frequency becomes 

u = k 2 + sk ul, (14) 

where u = Uq denotes the constant amplitude. Thus, with the plane wave ansatz, the effects 
due to the thermal correction (jSJ) vanish. 

Coherent modulational instabilities can be analyzed by letting u = {u(t) + 
Ui(t, z)) exp[ik z — iu (t)t] + c.c, where \u±\ <C u, the time variation of u is assumed slow 
compared to the time scale of the perturbation, and c.c. denotes the complex conjugate. 
Linearizing Eq. with respect to u±, dividing u\ into its real and imaginary part, and as- 
suming the wavenumber and frequency K and Q of the perturbation we find the dispersion 
relation 

= ~iDK 2 +2k K+2sKu 2 ± [s 2 K 2 u A + K 4 + 2sk K 2 u 2 - ADk K(Dk K + iK 2 + isk u 2 )] 1 

(15a) 

in the collisional case (j2J, and 

Q = -\iCKu 2 +2k K+2sKu 2 ± [s 2 K 2 u 4 + K 4 + 2sk Q K 2 u 2 - \CK 2 u 2 {\Cu 2 + ik Q + isu 2 )] 

(15b) 

in the collisionless case (JHJ). We note from ()15aj) and (J 15b)) that for s = 1 there is no 
instability and the perturbation modes are always damped. However, for s = — 1 this is not 
the case. Letting 7 = Im(f2) where 7 is the growth rate, we find from ()15a)) 

7 = _ DK i + (2k,K 2 u 2 - K 2 u 4 - if 4 ) 1/2 , (16) 

if we assume that D <C 1 and linearize for D. Similarly, in the collisionless case, the 
parameter C<1, and we may expand ()15b|) . keeping terms linear in C. Then we find 

7 = -\CKu 2 + (2k K 2 u 2 - K 2 u 4 - K A ) 1/2 . (17) 

Thus, the expressions (fTo]) and (|T7]) show that for s = — 1 4^ V s > Va we have growing 
modes due to the nonlinear evolution of the Alfven waves. However, these modes suffer 
damping due to dissipative processes. 

Next we will compare the above results for the coherent modulational instability with 
the case of partial coherent Alfven waves. 
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IV. EFFECTS OF PARTIAL COHERENCE 



Returning to Eq. we first look for solutions of the form p = p(t, k). In the collisional 
case (I12a|) Eq. (fTTj) takes the form d t p + 2Dk 2 p = 0. This is integrated to yield 

p = p (k)exp(-2Dk 2 t). (18) 

Equation ()18|) is in agreement with the monochromatic result (jl3b|) based on Eq. if 
we let po(k) = Io5(k — ko). In the collisionless case (j3J), the ansatz for p gives H = 0, 
where I(t) = J dkp(t,k). Thus, the expression (j!2b|) gives p = po(k) through Eq. (jllj) . in 
agreement with the monochromatic result. 

Next we look at perturbation around the background solutions. Thus, we let p = p(t, k) + 
Pi(k) exp(iKz — iflt) + c.c, where |pi| <C p an d \dtp\ ^ ^P, and linearize with respect 
to pi. With / = I(t) + hexp^Kz - iQt) + c.c. it follows that H w H + H 1 = H x = 
mil exp(iKz — iflt) + c.c. Thus, using 



and 



2 sin ( -Kd k ) p(t, k) = i [p(t, k + K/2) - p(t, k - K/2)} (19a) 



2 cos ( % -Kd k ) pit, k) = p{t, k + K/2) + p(t, k - K/2), (19b) 



2 

we find that the nonlinear dispersion relation 



[s-ia] f (k - K/2) pjt, k + K/2) - (k + K/2)pjt, k - K/2) 



2K J k + [sKI(t)-n-iL(k)]/2K 

valid for partially coherent Alfven waves. Here L(k) = DiK 2 /2 + 2k 2 ) and a = in 
the collisional case (J2J), and L{k) = and a = C/4 in the collisionless case 0. In the 
monochromatic case, we have po{k) = IoSik — ko). The dispersion relation f!2()l) then becomes 



Q = -\i(L+ + L_) + 2k K+ (2s - ia)KI (21) 
±{{s- ia) 2 K 2 I 2 + K A + 2{s - ta)k K 2 I - (L+ - L_) - L_) + ^ 2 + i(s - ia)k T\ } 1/2 

where L± = L(ko ± K/2). Thus, inserting the definitions of L{k) and a we regain the 
dispersion relations ()15|) when neglecting the slow background time variations. 

However, if the Alfven waves suffer random perturbations in e.g. the phase, the correlation 
between the waves may be nontrivial [14j. The partial coherence introduced by a random 
phase may be modeled by the Lorentzian distribution 1^] 



where A denotes the width of the distribution around ko- Due to the finite width of the 
Lorentz distribution the dispersion integral (j2*Uj) will have poles. In the collisional case (J2J) 
we obtain a rather lengthy analytical expression, where a complicated interplay between the 
spectral broadening, represented by A, and the collision parameter D takes place. Thus, we 
do not explicitly state this particular result here, but instead give the numerical solution of 
(120}) in conjunction with the Lorentzian distribution (}2*2*j) in the collisional case (J2J) (see Fig. 
[Q. However, in the collisionless case Q, the result is somewhat more compact and takes 
the form 



where f2 = — 2k K — sKI . 

Letting Q = Re(f2) + iT, where Re denotes the real part, we can solve the nonlinear 
dispersion relations numerically for the growth rate T as a function of the wavenumber K. 
The results are shown in Figs. Q]and|21 where we have used Iq = 0.5 and k = 1. We see that 
the parameter region with s = 1 exhibits new properties as compared to the coherent case. 
A new type of instability occurs due to the spectral broadening of the dispersive Alfven 
pump wave. 

V. SUMMARY 

We have investigated the nonlinear stability of broadband dispersive Alfven waves in 
magnetoplasmas, using a Wigner-Moyal formalism. A new mode structure with new insta- 
bilities, due to the finite spectral width of the dispersive Alfven pump wave, is found. The 
spectral broadening in conjunction with the kinetic modification of the DNLS equation can 
thus give rise to growing modes not present in the coherent case. 




-iKA 



AK 2 A 2 + K 2 ) 2 /2 + + K 2 f/2 + 2(Q + 2k K)(n - sKI ) 
[4A^ 2 A 2 + (fi - K 2 ) 2 ][4K 2 A 2 + + K 2 ) 2 } 



(23) 
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FIG. 1: The growth rate for different parameter values in the collisional case l|5j). The upper left 
panel has s = —1, while D = A = 0; the middle upper panel has s = —1 and D = 0, while A = 0.1; 
the right upper panel has s = — 1 and A = 0, while D = 0.1. The left lower panel has s = — 1 and 
D = A = 0.1, while the right panel has s = 1 and D = A = 0.1. Note that the last of these panels 
has no coherent counterpart. 
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FIG. 2: The growth rate for different parameter values in the collisional case ©. The upper left 
panel has s = —1, while C = A = 0; the middle upper panel has s = —1 and C = 0, while A = 0.1; 
the right upper panel has s = — 1 and A = 0, while C = 0.1. The left lower panel has s = — 1 and 
C = A = 0.1, while the right panel has s = 1 and C = A = 0.1. Note that the last of these panels 
has no coherent counterpart. 
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